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We present a model for weak CC induced nuclear reactions at energies of interest for current and 
future neutrino oscillation experiments. This model is a natural extension of the work of Refs. 
where the QE contribution to the inclusive electron and neutrino scattering on nuclei was analyzed. 
The model is based on a systematic many body expansion of the gauge boson absorption modes 
that includes one, two and even three body mechanisms, as well as the excitation of A isobars. The 
' whole scheme has no free parameters, besides those previously adjusted to the weak pion production 

^sQ I off the nucleon cross sections in the deuteron, since all nuclear effects were set up in previous studies 

, of photon, electron and pion interactions with nuclei. We have discussed at length the recent 

'"JT ' charged current quasi-elastic MiniBooNE cross section data, and showed that two nucleon knockout 
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'"17. I. INTRODUCTION 

Ph. 

' ' The interaction of neutrinos with nuclei at intermediate energies plays an important role in the precise determination 
of neutrino properties such as their masses and mixing parameters. It can also provide relevant information on the axial 
i hadronic currents. The statistical significance of the experiments is rapidly improving. However, the data analysis 
needs to consider a large number of nuclear effects that distort the signals and produce new sources of background 
that are absent in the elementary neutrino nucleon processes. 
^ In this context, it is clearly of interest the elaboration of a theoretically well founded and unified framework in 

J — i which the electroweak interactions with nuclei could be systematically studied. Furthermore, the recent measurements 
f — ' of the cross sections for several channels provide a serious benchmark to the theoretical models. An excellent 
f — ■ ^ review of the current situation can be found in Ref . [TJ . 

CN ■ A suitable theoretical model should include, at least, three kinds of contributions: (i) quasielastic (QE) for low energy 
I transfers, (ii) pion production and two-body processes from the QE region to that around the A(1232) resonance peak, 
. and (iii) double pion production and higher nucleon resonance degrees of freedom induced processes at even higher 
T-H ■ energies. A word of caution is needed here, because the same words could refer to somehow different magnitudes in 
, the literature. For instance, whereas in most theoretical works QE is used for processes where the gauge boson or 
^ ■ Z° is absorbed by just one nucleon, which together with a lepton is emitted^, in the recent MiniBooNE papers, QE is 
related to processes in which only a muon is detected. This latter definition could make sense because ejected nucleons 
are not detected in that experiment, but includes multinucleon processes and others like pion production followed by 
^ I absorption. However, it discards pions coming off the nucleus, since they will give rise to additional leptons after their 
■ " " decay. In any case, their experimental results cannot be directly compared to most previous calculations. 

The QE processes have been abundantly studied. Simple approaches using a global Fermi gas for the nucleons 
and the impulse approximation are good enough to describe qualitatively electron scattering but more sophisticated 
treatments of the nuclear effects are necessary to get a detailed agreement with data. There are different kinds of 
models like those based on the use of proper nucleon spectral functions 8 10], others in which nucleons are treated in 



■ Il2l a nd models based on a local Fermi gas including many body effects such as spectral 
Concerning the elementary process, -\- N ^ I -\- N', the hadronic vector current 



a relativistic mean field [11 
functions [l^ and RPA 

is well known from electron scattering. The axial current, after the use of the partial conservation of the axial current 
to relate the two form factors and assuming a dipole form, depends on two parameters: gA, that can be fixed from 
the neutron /3 decay and the axial mass Ma- The value of Ma established from QE data on deuterium targets is 



^ This follows the traditional nomenclature of electronuclear scattering. Note that in same cases the resulting nucleon after the absorption 
of the gauge boson is not emitted, but rather it could be trapped and form part of a bound state of the daughter nucleus (discrete 
transition). 
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Ma = 1.016 ± 0.026 [13] GeV. A consistent result is obtained from n electro-production after chiral corrections are 
incorporated [Tsl. [Toj. 

The predicted cross sections for QE scattering are very similar for most models. See, e.g., the compilation shown 
in Fig. 2 of Ref. ^2d\ . On the other hand, the theoretical results are clearly below the recently published MiniBooNE 
data The discrepancy is large enough to provoke much debate and theoretical attention. Some works try to 
understand these new data in terms of a larger value of Ma- For instance, in Ref. B a value of AIa — 1.35 ± 0.17, 
that also fits the shape, is suggested. Consistent values are obtained in Refs. |2l|- [23 . This idea is not only difficult 
to understand theoretically, but is also in confiict with higher energy NOMAD data 12411 {Ma = 1.06 ± 0.02[stat) ± 
0.06{syst) GeV). In another line of research, the role of meson exchange currents 12511 and superscaling [2g| have 
been also estimated recently. Finally, another idea has been explored in Refs. [13, l28j |. which include two nucleon 
mechanisms (and others related to A excitation) and reproduce MiniBooNE QE data without the need of a large value 
of Ma- These latter results suggest that much of the experimental cross section can be attributed to processes that 
are not properly QE, stressing again the need of a unified framework dealing with all relevant mechanisms, namely tt 
production and multinucleon excitation. 

The matter of tt production induced by neutrinos is also of much interest ^29-3^]- The elementary reaction on 
the nucleon, at low and intermediate energies, includes both background and resonant mechanisms. The background 
terms can be obtained from the chiral lagrangians. The resonant terms contain some free parameters that have 
been adjusted to ANL and/or BNL old bubble chamber data. Still, the experimental data have large normalization 
uncertainties which are certainly reflected in the theoretical models. At low energies, the A(1232) resonance plays 
a very important role in this process, and for small values only one form factor (C^) is relevant. Thus, special 
attention has been addressed to its study with recent results ranging from C^(0) = 1.19±0.08 [3l|, obtained neglecting 
the non resonant background, to C^(0) = 1.00 ± 0.11 [sH] in a more complete model. This latter value implies a 20% 
reduction with respect to the off-diagonal Goldberger-Treiman relation (GTR). In nuclei, several effects are expected 
to be important for the tt production reaction. First, the elementary process is modified by Fermi motion, by Pauli 
blocking and more importantly by the changes of the spectral function of the A resonance in the medium. In addition, 
the final pion can be absorbed or scattered by one or more nucleons. This latter kind of effects do not modify the 
inclusive neutrino nucleus cross section and thus are out of the scope of this paper. 

Our aim in this work is to extend the model of Ref. , which studied QE scattering. We will include two nucleon 
processes and tt p roduction in a well established framework that has been tested, for instance, in electron and photon 
scattering [l|, l37l |. This will extend the range of applicability of the model to higher transferred energies (and thus 
higher neutrino energies) and allow for the comparison with inclusive data which include the QE peak, the A resonance 
peak and also the dip region between them. The structure of the paper is as follows: In Sect.|lTl we start establishing 
the formalism and reviewing briefly the approach for QE scattering of Ref. [^. Then, we consider pion production 
mechanisms and two nucleon processes. Next, we discuss with special care the role of the A resonance and how it is 
affected by the nuclear medium. In Sect. |in]we present and discuss some of the results derived from the model, and 
compare these to the recent MiniBooNE charged current (CC) QE and SciBooNE total cross section data. Finally in 
Sect. IIVI we draw the main conclusions of this work. 



II. CC NEUTRINO/ANTINEUTRINO INCLUSIVE NUCLEAR REACTIONS 

A. General considerations 



We will focus on the inclusive nuclear reaction 

i^iik) + Az ^ r ik') + X (1) 

driven by the electroweak CC. The generalization of the obtained expressions to antineutrino induced reactions 
is straightforward (see Subsect. IIIFI ). The double differential cross section, with respect to the outgoing lepton 
kinematical variables, for the process of Eq. ([1]) is given in the Laboratory (LAB) frame by 

. M^L„H- (2) 

dVL{k')dE[ |fc I 47r2 

with k and k' the LAB lepton momenta, E[ = (P^ + m^)^/^ and m; the energy and the mass of the outgoing lepton, 
G = 1.1664 X 10^^^ MeV^2, the Fermi constant and L and W the leptonic and hadronic tensors, respectively. To 
obtain Eq. Q we have neglected the four-momentum carried out by the intermediate VF— boson with respect to its 
mass, and have used the relation between the gauge weak coupling constant, g = e/ sinOw, and the Fermi constant: 
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= /SM'^r, with e the electron charge, 6w the Weinberg angle and Mw the W— boson mass. The leptonic 
tensor is given by^: 

Lfia = -^pcr + i-^pcr ~ ^^^a ^cr^p ~ 9fJ.a-k ■ k + ie^crQ^fc . (3) 

The hadronic tensor corresponds to the charged electroweak transitions of the target nucleus, i, to all possible final 
states. It is given by^ 

with P the four-momentum of the initial nucleus. Mi = V the target nucleus mass, P'^ the total four momentum 
of the hadronic state / and q — k ~ k' the four momentum transferred to the nucleus. The bar over the sum denotes 
the average over initial spins. The hadronic CC is given by 

j^, = ^uT^ll - 75)(cos0c*d + sin 0c*.) (5) 

with 'ifd and ^'s quark fields, and 9c the Cabibbo angle. By construction, the hadronic tensor accomplishes 

W^"" = Wr + iW^" (6) 

with W^'^ i^a'^) I'eal symmetric (antisymmetric) tensors. The hadronic tensor is determined by the boson 
selfenergy, 11^ (g), in the nuclear medium. We follow here the formalism of Ref. where it is shown that within 
the local density approximation, the hadronic tensor can be written as 



Wr = -e(g°) r-^] / ^ Im [n^:,- + (g; p) (7) 




Wr = -0(9°) ( ^ 1 / — Re [n;V" - n^^] (g; p). (8) 



Then, the differential cross section for the reaction in Eq. ([T]) is given by 
d'a^i \k'\ / 2V2 \ fd^r 




dn{k')dk'0 \k \ 2Tr^ y g J J 2Tr 

\k'\ ( 2^/2V fd^r 



Im 



L^,n;!^(g;p) e(g") 



with 8(...) the Heaviside step function. 

The in medium gauge boson (M^*) selfenergy depends on the nuclear density p{r). We propose a many body 
expansion for 11^, where the relevant gauge boson absorption modes would be systematically incorporated: absorption 
by one nucleon, or a pair of nucleons or even three nucleon mechanisms, real and virtual meson (tt, p, • • •) production, 
excitation of A of higher resonance degrees of freedom, etc. In addition, nuclear effects such as RPA or Short Range 
Correlations (SRC) will be also taken into account. Some of the basic W^— absorption modes are depicted in Fig. [TJ 



B. Quasielastic scattering 

The virtual W'^ can be absorbed by one nucleon leading to the QE contribution of the nuclear response function. 
Such a process corresponds to a one particle-one hole (Iplh) nuclear excitation (first of the diagrams depicted in 
Fig. ID). This contribution was studied in detail in Ref. |2)|^. Here, we will just briefly discuss the main features of the 
model. In Ref. starting from a Local Fermi Gas (LEG) picture of the nucleus, which automatically accounts for 
Pauli blocking and Fermi motion, several nuclear corrections were incorporated, among others: 



^ We take eoi23 = -1-1 and the metric g^" = (+, — , — , — ). 

^ In Eq. ^ the states are normahzed such that {p\p') = {2n)'^2poS^ {p — p') and the sum over final states / includes an integration 

/d^ p ■ 
(27r)^2S ' ^^"-^ particle j making up the system /, as well as a sum over all spins involved. 

^ The extension of this scheme for Neutral Currents (NC) was discussed in Ref. \Ti . 
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FIG. 1: Diagrammatic representation of some mechanisms contributing to the W+— selfenergy. 



• A correct energy balance, using the experimental values, was enforced. 

• Coulomb distortion of the charged leptons, important at low energies, was implemented by using the so called 
"modified effective momentum approximation" . 

• Medium polarization (RPA), including A— hole degrees of freedom and explicit pion and rho exchanges in the 
vector-isovector channel of the effective nucleon-nucleon force, and SRC effects were computed. 

• The nucleon propagators were dressed in the nuclear medium, which amounts to work with nucleon spectral 
functions (a LFG of interacting nucleons) and it also accounts for some reaction mechanisms where the gauge 
boson is absorbed by two nucleons. 



This model is a natural extension of previous studies on electron ,1J, photon ^37j and pion |38l - l42| dynamics in nuclei. 
Even though the scarce existing CC data involve very low nuclear excitation energies, for which specific details of 
the nuclear structure might play an important role, the model of Ref. provides one of the best existing combined 
descriptions of the inclusive muon capture in '^^C and of the ^^C (i/^, n~)X and ^^C (i^e, e~)X reactions near threshold. 
Inclusive muon capture from other nuclei is also successfully described. 

The theoretical errors affecting the predictions of Ref. [1] were discussed in Ref. [i^. There, it is concluded that is 
sound to assume errors of about 10-15% on the QE neutrino-nucleus (differential and integrated) cross section results 
of Ref. [2^. 

The LFG description of the nucleus has been shown to be well suited for inclusive processes and nuclear excitation 
energies of around 100 MeV or higher [l], ISTl - li^ . The reason is that in these circumstances one should sum up over 
several nuclear configurations, both in the discrete and in the continuum. This inclusive sum is almost insensitive to 
the details of the nuclear wave function^ , in sharp contrast to what happens in the case of exclusive processes where 



^ The results of Ref. Q| for the inclusive muon capture in nuclei through the whole periodic table, where the capture widths vary from 
about 4x10* s-l in to 1300 xlO"* s"! in '^°^Ph, and of the LSND measurements of the (j/^,/^-)^ and (ue,e-)X reactions 
near threshold indicate that the predictions of our scheme, for totally integrated inclusive observables, could be extended to much 
smaller, 10 or 20 MeV), nuclear excitation energies. In this respect, Refs. [43 | and [45l | for inclusive muon capture and radiative 
pion capture in nuclei, respectively, are enlightening. In these works, continuum shell model and LFG model results are compared for 
several nuclei from ^-^C to ^"^Pb. The differential decay width shapes predicted for the two models are substantially different. Shell 
model distributions present discrete contributions and in the continuum appear sharp scattering resonances. Despite the fact that those 
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FIG. 2: Left: Definition of the kinematical variables used in tfiis work. Rigiit: Model for the W^N N'tt reaction. It 
consists of seven diagrams: Direct and crossed A(1232)— (first row) and nucleon (second row) pole terms, contact and pion 
pole contribution (third row) and finally the pion-in-fiight term. Throughout this work, we will label these contributions by: 
AP, CAP, NP, CNP, CT, PP and PF, respectively. The circle in the diagrams stands for the weak transition vertex. 



the final nucleus is left in a determined nuclear level. On the other hand, the LFG description of the nucleus allows 
for an accurate treatment of the dynamics of the elementary processes (interaction of gauge bosons with nucleons, 
nucleon resonances, and mesons, interaction between nucleons or between mesons and nucleons, etc.) which occur 
inside the nuclear medium. Within a finite nuclei scenario, such a treatment becomes hard to implement, and often 
the dynamics is simplified in order to deal with more elaborated nuclear wave functions. 



C. The virtual self-energy in pion production: Iplhlyr contribution 

In this subsection, we calculate the contribution to the cross section from W'^ gauge boson self-energy diagrams 
which contains pion production in the intermediate states. We will use the model for the CC neutrino-pion production 
reaction off the nucleon, 

Pi{k) + N{p) ^ r (fc') + Nip') + ^{K) (10) 

derived in Refs. [H, [s^]. This process, at intermediate energies, is traditionally described in the literature by means 
of the weak excitation of the A(1232) resonance and its subsequent decay into Nit. In Ref. [11], some background 
terms required by the pattern of spontaneous chiral symmetry breaking of QCD are also included. Their contributions 
are sizable and lead to significant effects in total and partially integrated pion production cross sections even at the 
A(1232)— resonance peak, and they are dominant near pion threshold. The model consists of seven diagrams (right 
panel of Fig. [2]). The contributions of the different diagrams are calculated by using the effective Lagrangian of the 
SU(2) nonlinear cr— model, supplemented with some form-factors (see Ref. [32j for details). In this work, we will use 
the set IV of form factors compiled in Table I of Ref. j35j] . The available data set on neutrino and antineutrino pion 
production on nucleons is described reasonably well. Nonetheless, we must mention, that the experimental data still 
have large uncertainties and there exist conflicting data for some channels. 

The discussed model can be considered an extension of that developed in Ref. [l| for the eiV — > e'N-K reaction. 
For the latter case, the model, that contains a theoretically well founded description of the background amplitudes, 
provides the same level of accuracy [S^ as the MAID model [i^ , which ensures its applicability to the leptoproduction 
processes at least up to < 1.4 GeV, being W the outgoing ttN invariant mass. 



distinctive features do not appear in the LFG differential decay widths, the totally integrated widths (inclusive observable) obtained 
from both descriptions of the process do not differ in more than 5 or 10%. The typical nuclear excitation energies in muon and radiative 
pion capture in nuclei are small, of the order of 20 MeV, and thus one expects that at higher excitation energies, where one should sum 
up over a larger number of nuclear final states, the LFG predictions for inclusive observables would become even more reliable. 
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FIG. 3: W^— self-energy obtained by folding the WN — >■ vrA'^' amplitude (A is the charge of the pion). 



We move now to the computation of the W'^ gauge boson self-energy diagrams which contain pion production in 
the intermediate states. This is readily accomplished by taking the W~^N — >■ irN' amplitude of Fig. [5] and folding it 
with itself. One gets then the diagram of Fig. [3] where the circle stands for any of the 7 terms of the elementary model 
for WN — > ttN'. The solid lines going up and down in Fig. [3] follow the standard many body nomenclature and stand 
for particle and hole states respectively. The W^— self-energy corresponding to this diagram (actually 49 diagrams) is 
readily evaluated and gives® (we will label it as IplhlTr) 

X Tr((;} + A/)70jf70(/ + M)jl) (11) 



E 



X 



X D 



2V2J {2^YJ {2^f2E{p)2E{p + q-K)q'^-kl + E{p)-E{p')+it 

.(fc.)Tr {{i> + M)7"jf 70(/ + M)j\) + [{q - fc.) ^ -(g - fc.)] (12) 

where p' ^ p + q — k-^, j'^ is the amputated amplitude^ for the N — > N'tt^ process, which is obtained by summing 
up the contributions of all diagrams of the right panel of Fig. The contribution to of each diagram is given by 
their relation to the full amplitudes given in Eq. (51) of Ref. [33 |. 



Jcc+ 



u{p')j'XiP,q,P =P + q-k^,k^)u{p), AP,CAP,NP,CNP,CT,PP,PE (13) 



The indices N,N' in Eq. stand for the hole and particle nucleon states respectively and UNip) = Q{kp — \p\) 
is the occupation number in the Fermi local sea, with kp{r) = (37r^pAr(r))^/^ and pn^i") the density of nucleons of a 
particular species N = n or p (p(r) — Pp{r) + p„(r)), normalized to the number of protons or neutrons. Besides, E(p) 
is the energy of the nucleon y/p'^ + A'P — kp/2M, with M its mass and kpir) — (37r^p(r)/2)^/'^, and Dyr is the pion 
propagator 



1 

kl-' 



DAk.) - 72 r— (14) 



with m^r the mass of the pion. Besides, the nucleon propagator in the medium reads, 

S{p;p) = {^ + M)G{p;p) (15) 



® In Eq. jlTJ, it is necessary to subtract the free space contribution, i.e., the one that survives for vanishing nuclear densities. This 
contribution will renormalize free space couplings and masses. To obtain Eq. 1121 1. we have neglected the contribution of the antiparticle 
pole (p" = —E{p) — it) in the integration, this also gets rid of the vacuum contribution that needed to be subtracted. 
The dependence of j'^ on the A'^, N' , A channel is understood and it is not made explicit. 
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Gip;p) 



A'P 



n{p) 



1 - n{p) 



E{p) + ie \p° - E{p) - ie pO - S(p ) + ie 



(16) 
(17) 



A further simplification can be done by evaluating the amplitudes at an average momentum, which allows to 
take the spin trace in Eq. ([T2|) out of the d^p integration. This latter integration can be now done, and it gives, up 

to some constants, the Lindhard function, U R{q — fc^, fc^, kp ) (see appendix B of Ref. Q)- We take {\p\) = \J\kp 
and a direction orthogonal to the plane defined by the pion and the virtual gauge boson. Within this approximation, 
we find 



inw;lphl7r(9) 



2^/2 / 4M2 



E 

N.N'X 



(2^ 



■D^{k^)UR{q - fc,, k^, kf)Af''' [{p),q,p' ={p)+q- Kfl^) 



A^" = iTr(((|}>+M)7''0f)7°((^^)+^-|^. + A/)01)) 



(19) 



where (j^) stands for calculated with the average hole momentum (p). To find the contribution to the hadron 
tensor W^'^ of the many body diagrams depicted in Fig. [3l we remind that by construction 



with A^'^ (^a*^) real symmetric (antisymmetric) tensors, and thus 



Im 



-^^W;lphl7r 



^^W;lphl7r 



Re 



^^W;lphl7r 



^^W;lphl7r 



(20) 



(21) 



where li'^.^^^^^ 



s(a) 



is defined as in Eg. p^ . but replacing the full tensor A^'^ by its symmetric (antisymmetric) A^J^ 



parts. The imaginary part of H'"^ 



VK;lphl7r 



s(a) 



(a) 

can be obtained by following the prescription of the Cutkosky's rules. In 



this case we cut with a straight horizontal line the intermediate particle and hole states and the pion. Those states are 
then placed on shell by taking the imaginary part of the propagator. Technically the rules to obtain Imll^.-^pj^j^^ 
reduce to making the substitutions: 



nlV^(g) 



.^y^^ ^ 2iImn^V"(Q)e(g°) 
D^{k^) -> 2\lu,D,{k^)Q{kl) = -2T:i5{kl - ml)Q{kl) 
UR{q-k^,k^,k^') ^ 2ilmUR{q-k^,k^,k^')Q{q° ^k"^) 



Thus, we readily obtain 



<phi.(9) = -e(g") 



2Af2 



^ E f-t ^^^\,-.''°^ lmUn{q-K,k^,k^')A-^^' 



(2^)3 2w{k^) 



(a) 



(22) 
(23) 
(24) 



(25) 



with iolk^r ) the pion on-shell energy. The approximation done saves a considerable amount of computational time 
since there are analytical expressions for lmUu{q — fc^, kp , kp ) 2]^. 



In the small density limit ImU H{q, , k^' ) ~ --rrp^MS [^q'^ + M - ^ AP + q'^j I \J M'^ + Substituting this into Eq. ([25j one 
obtains 

lim W^",, ~ + + + l^w+ +1 (26) 

being Z and A'^ the number of protons and neutrons of the nucleus, and ^'^j^ ^j^, \ the hadronic tensor for CC pion production off 
the nucleon (see Eq. (4) of Ref. [32ll 'l. In this way, the strict impulse approximation is recovered. By performing the integral in Eq. II25I I. 
one accounts for Pauli blocking and for Fermi motion. 
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FIG. 4: VP'— self-energy obtained from the one in Fig. |3]when the pion line is allowed to excite a particle-hole. 



D. The dip region: 2p2h absorption 

In the two previous subsections we have discussed the dominant contributions to the inclusive v^v nucleus CC cross 
section at low energies, namely QE scattering and pion production. In this subsection, we present a model for 2p2h 
mechanisms, which could be very relevant in the description of the region of transferred energies above the quasielastic 
and below the A— resonance peaks (the dip region). 



1. 2p2h mechanisms driven by the longitudinal part of the effective spin-isospin ph-ph interaction 

Let us consider again to the generic diagram of weak pion production of Fig. [3] and allow the pion to excite a 
particle-hole. This leads us to the diagram of Fig. |4l This is still a generic diagram which actually contains 49 
diagrams when in the shaded circle we put each one of the terms of the WN — ^ ttN amplitude of Fig. [5] The 
diagrams in Fig. 2] contribute to Imll^^ according to Cutkosky rules when they are cut by a horizontal line, and the 
2p2h are placed on shell. The contribution of the diagram of Fig. |4] can be obtained from that of Fig. [3l given in 
Eq. ([HI), by replacing (A^ = 1.2 GeV, J^nn/'^^ = 0.08) 

f2 ^ A2 _ 7772 

D^ik^) ^ Dlik,)F^{k^)l^klU,{k,), F^{k) = I (27) 

where Ux is the Lindhard function for a particle-hole excitation by an object of charge A: this is, twice or C/]^'" 
for the excitation by a charged pion or f/^'^ + C/^'" for the excitation by a neutral pion. The pion form factor F^{kT^) 
appears because now the pions are off shell. 

We can again simplify the expression by taking an average nucleon momentum of the Fermi sea to evaluate the 
amputated amplitudes for the W^N — >■ N'tt^ process. This allows us to factorize the Lindhard function and following 
the prescription of the Cutkosky's rules we get, 

W,';2^il) = ^(1°)^ [ ^e(<z°-fc°)F2(fc.)ImC7^(q-fc,,fc^,fcr)A''^x 

xDlik^) Fi (fc, ) klQ{kl)lmUx{k^) (28) 

Next, we have implemented several improvements that account for well established many body corrections: 
1. In the above expression of Eq. (pSj) . we have replaced 

Dl{k^)F^,{k^)^-^Lkli^Ux{k.) ^ Im {- / - -^^ ^T^u(^ M2 (29) 

\k4 — rn^ — II(fcTr)/ \ki — mi — Vl[kT^)Y 



where for the selfenergy of a pion of charge A, we have taken [4l| 



n(fc.) = Fl{k^)lkfkl ^ (30) 



FIG. 5: W^— selfenergy diagram obtained from the QE Iplh excitation term (first of the diagrams depicted in Fig.[T]) by dressing 
up the nucleon propagator of the particle state in the ph excitation. 



where 

U{K)^UN{K) + UA{k^) (31) 

is the non-relativistic Lindhard function for ph + Ah excitations^ including direct and crossed bubbles [U 13 ' '^^ 
contrast to Ur which only contains the direct bubble of a particle-hole excitation (the only one which contributes 
to Im[/jv for > 0). When evaluating Inin(fc7r) in the numerator of Eq. (f29|) we have not considered the part 
that arises from putting the Ah excitation on-shell that would correspond to a 2p2h+l7r mechanism. We 
expect this latter contribution to be small at the considered energies. Note also, that by using U to compute 
the pion selfenergy, we have neglected small relativistic and pp ^ p„ corrections. By means of Eq. (p9)) . we 
have implemented the Dyson re-summation of the pion selfenergy, and have improved on this latter one by 
incorporating the Lorentz-Lorenz effect, driven by the short range Landau Migdal parameter g' (40l |. going 
in this way beyond Iplh excitation^" in the evaluation of n(fc7r). We have used g' = 0.63, as in previous 
works [il,i,|43,|4ll. 

2. Let us now pay attention to the diagram of Fig. [SJ which is already implicit in the generic diagram of Fig. 2] 
when the NP amputated amplitude is considered in both weak vertices. This VF— selfenergy contribution can 
be obtained from the QE Iplh excitation term (first of the diagrams depicted in Fig. [1} by dressing up the 
nucleon propagator of the particle state in the ph excitation. Indeed, this, among other contributions, was 
already taken into account in the QE study carried out in our pre vious work of Ref. since there dressed 
nucleon propagators deduced from a realistic nucleon selfenergy [49| were used. To avoid double counting, we 
subtract the contribution of the NP-NP diagram of Fig. [5] from Eq. 



3. When in one of the weak vertices of Fig. HI the NP term is considered, the prescription of taking an average 
nucleon momentum of the Fermi sea used to obtain Eq. (l28l) turns out to be not appropriated. The reason is 
that when placing the 2p2h excitation on shell, through Cutkosky rules, we still have the nucleon propagator 

with momentum p + q (this is part of the amputated amplitude ) . This propagator can be still placed 

on shell for a virtual W and thus, there exists a single pole in the (Pp integration^^. In this situation, one can 
not take an average for p, as we have implicitly assumed in Eq. (j28p . and we have improved such prescription 
as follows. In this latter equation, it appeared the tensor A^f^, which in turn is defined in Eq. (|T9| by using an 
average for the hole three momentum pto calculate both, the amputated WN — >■ Ntt amplitudes and ^, that 
also appears in the trace that defines A'^'^. Instead of this, we have computed an average of the whole trace. To 
this end, we have numerically performed the integral over the angle formed by p and q, using still an average for 
the modulus of p and taking this momentum in the XZ plane (recall that q defines the Z-axis) . All pathologies 



^ The functions Um and Ua are defined, e.g., in Eqs. (2.9) and (3.4) of Ref. [4^ . Un incorporates a factor two of isospin with respect to 
Uji, such that Im[/jv = SlmC/^ for symmetric nuclear matter, up to relativistic corrections. 

It corresponds to replace the ph excitation of the right-hand in Fig. |4]by a series of RPA excitations through ph and Ah excitations, 
driven by the longitudinal part of the effective spin-isospin interaction. In Subsect. Ill D 21 we do something similar for the case of 
2p2h mechanisms driven by p— meson exchange, and there we show graphically in Fig. llOl the RPA series, in that case induced by the 
transverse part of the effective spin-isospin interaction. 

This cut will also contribute to the nuclear response to the weak probe. But, while it will aff'ect to the QE region, it is expected to 
be small and considerably difficult to calculate from the computational point of view (see Eq. (80) of Ref. [Ij). Thus, for the sake of 
simplicity we have not considered such contribution. 
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arise from the p + q nucleon propagator hidden in the amputated amplitudes, which can be put on the mass 
shell, and thus the contribution of these diagrams depends critically on the angle formed by p and q, while it 
shows a very smooth dependence on the rest of kinematical variables of the hole momentum p. Thanks to the 
approximations of using an average for the modulus of p and fixing the {p, (f)— plane , we avoid to perform two 
nested integrals, with the obvious benefit in computation time. We have checked that the results are accurate 
at the level of 5-10%. To be more specific, in Eq. ([SS)) . we have replaced A'^'-'- by 



(32) 



with fi = q- p/\q\\p\. To speed up the numerical integration, we have also given a small width (^ 10 MeV) to 
the p + q nucleon. Results do not depend significantly on this choice. 

For consistency, we have also performed this angular average for all contributions implicit in Fig. |4l though 
the prescription of using an average for p leads to accurate results in all cases except those involving the NP 
amputated amplitude discussed above. 

4. In the terms involving the NP, amputated amplitude (interferences with the rest of amplitudes of Fig.[T]), there 
always appears a pion emitted after the WN vertex that couples to the second ph excitation (see for instance the 
line labeled as tt in Fig. There, one is assuming a pion exchange interaction among the two ph excitations. 
We have improved on that, and have replaced it by an effective longitudinal interaction, Vi, 

Viik) = l^LlF^{k)-^+gl{k)], (33) 




which besides pion exchange includes SRC driven by the Landau Migdal parameter g'i{k) (see Refs. |4Cl I4lll48l|). 
To achieve this, we have multiplied the amputated amplitude jj^^^ by a suitable factor, 

Ja =^3 a X I 1 + ^ I (34) 

We have taken the same prescription also for those terms that include the CNP, AP and CAP amputated 
amplitudes. 

We have also considered the transverse channel interaction, Vt, 

Vt{k) = 

of the effective spin-isospin interaction among the two ph excitations. Here, nip = 0.77 GeV. The SRC 
functions g'l and g[ have a smooth fc— dependence [S^H^, which we will not consider here^^^ and thus we will 
take g[{k) = g'i{k) — g' — 0.63, as it was done in the study of inclusive nuclear electron scattering carried out in 
Ref. [l| , and also in the previous work on the QE region of Ref. . To account for such contribution to the 2p2h 
absorption cross section is slightly more complicated, because the tensor A'^'^ does not account for p— meson 
production in the primarily weak vertex. Details will be given in Subsect. HID 21 

The cut which places the two ph on shell in the diagrams of Fig. |4]is not the only possible one. In Fig.[7l we show 
a different cut (dotted line) which places one ph and the pion on shell. As done for real 37] and virtual [1] photons, 
we neglect this contribution in the non resonant terms, because at low energies where these pieces are important, the 
(W, tt) channel is small and at high energies where the {W, tt) contribution is important, this channel is dominated by 
the A excitation and there this correction will be properly incorporated. 

We have also considered two body diagrams, where each W boson couples to different bubbles (Fig. [5]). Its 
contribution to the hadron tensor, taking average momenta for both hole nucleon momenta in first place, reads 



This is justified because taking into account the fe— dependence leads to minor changes for low and intermediate energies and momenta, 
where this efi'ective ph-ph interaction should be used. 



11 




FIG. 6: W^— selfenergy diagrams in which one of the vertices contains the NP term of the WN Nn amphtude, while the 
other one (filled circle) contains all terms except that one. 




FIG. 7: Same as Fig. U showing the cut which places one particle-hole and the pion on shell. 



'^'^2p2h-2b 



^2b 



2V2M 



1 



-Tr 



+ M) i^:^ - i) 75 + M) {f^ 



1. 



-Tr 



{W+n pn°) {W+n mT+) - {W+n n7r+) {W+n pn°) 



(36) 



+ {W+n -> p7r°) {W+p 



(37) 



To compute this moderately small term, we have taken a proton-neutron symmetric Fermi sea, i.e., = kp. We have 
labeled by p and I, the left and right bubble hole four momenta, and to compute (j^) and (i^^), the pion momenta 
are fc^r and /ctt ~ q, respectively. Besides, we write explicitly in Eq. (1371) the sum over all charge combinations. 




FIG. 8: Two particle-two hole W— selfenergy Feynman diagram where the outgoing gauge boson couples to the second nucleon. 
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Finally, we have improved on our evaluation of Vr^''2h--2b given in Eq. ((37|) by implementing also here the refinements 
3 and 4, previously discussed, devoted to the improvements in the computation of W2p2h- However, as in Ref. [1] for 
the inclusive electron-nucleus case, we have only considered the contributions stemming from the longitudinal part 
of the effective spin-isospin ph-ph interaction, driven by pion exchange (depicted in Fig. [S]) , and we have neglected 
those induced by the transverse one. 

2. 2p2h mechanisms dnven by the transverse part of the effective spin-isospin ph-ph interaction 

In this subsection we evaluate the contribution to the weak nuclear response of the 2p2h absorption terms driven by 
the transversepart of the effective spin-isospin ph-ph potential used in previous studies on electron [ij, photon (3^ 
and pion [38l - l43 | interaction with nuclei. In the model of Ref. [s!], this transverse interaction arises from p— exchange 
modulated by SRC. The major difficulty here, as compared with the previous works mentioned above, arises from 
the fact that we are using a relativistic description of the weak transition process. Thus, the first step is to model 
NNp and NAp relativistic Lagrangians, which give rise, in the non-relativistic limit, to the transverse potential of 
Eq. (jSSp . A convenient set of interaction Lagrangians is, 

CnAp = -i^^v^4'.757p7'^*(5^r-5V")+h.c. (38) 

where ^ = ^ ^ J is the nucleon field, p'^ is the p— meson Proca field^^, is a Rarita Schwinger J"" = 3/2+ field, ft 

is the isospin transition operator^'*, f are the isospin Pauli matrices and /*jva ~ 2.14. 

Next, we consider the W^N — > N' p process and find the NP, CNP, AP and CAP amputated amplitudes obtained 
from the above Lagrangians. The Feynman diagrams for these four amplitudes are like those depicted in the right 
panel of Fig. [2]by replacing the outgoing pion by a p— meson. We will denote the amputated amplitudes by t'^. They 
are defined by 



''CC+ 



^u{p'K"{p,q,p' ^P + q~kp,kp)u{p)eUkp), t = AP,CAP, NP,CNP (39) 
with e* , the p— meson polarization vector and 



finds that t'^^^ 



can be obtained from the W'^N -> N'tt amplitudes j^^^ 



the full W^N — !• N'p amplitude for each mechanism. One readily 
, given in Eq. (51) of Ref. 32], with the 



following replacements deduced from the appropriate meson- A^A^^ and meson- A^A vertices, 

NP and CNP terms : |^^75 ^ V^/c'V^^e*'" 

AP term : k^Pc.p#" ^ i ^C^Ts {^Pa0 - k^laP^p) ¥e*^ 
CAP term: #"fcfP„0 ^ -iv^#" (Pa.75^ - i'ar,A:"757.) e*" (40) 



where k is now the p— meson momentum {p + q = p' + k). Now, we can compute the contribution of the 16 diagrams of 
Fig.[9]to the inclusive neutrino-nucleus cross section, when the two ph excitations are put on shell. The corresponding 
contribution to the hadron tensor reads, with the approximations discussed in the previous sections in order to 
factorize two Lindhard functions, 

nV*!?) = -^il")^2 J ^,e{q'~k')F^{k)lmUn{q~k,k^,kf )B'''^x 



f2 



xDlik)F^ik)Cll^kle{k'')lmUxik) (41) 



Here = (p^" — ip2)/v^ creates a p from the vacuum or annihilates a p+ and the field creates or annihilates a p'^ . 
It is a vector under isospin rotations and its Wigner-Eckart irreducible matrix element is taken to be one. 
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FIG. 9: Two particle-two hole Vl^— selfenergy Feynman diagram driven by p— exchange. The cut (dotted line) that places the 
2p2h on-shell is also displayed. The empty circle contains the direct and crossed nucleon and A— pole terms of the WN Np 
amplitude. 

with^'^ Dp{k) ~ l/(fc^ ~ "ip)- The form factor Fp{k), Eq. ([55]) . appears because the p's are off shell. Also here, when 
placing the 2p2h excitations on shell, we have that the nucleon propagator with momentum p + q (this is part of the 

amputated amplitude t'^ ) can be placed on shell for a virtual W. Thus, as discussed above in Eq. ([5^ . we define 

the tensor B'^^ as an angular average of the traces that appear in the evaluation of the diagram. Namely, 

5-^^ = 1 l^\f,^Tr{{^ + M)j°{tArj7''i^ + 4-f^ + M){tAr"} (42) 

with fi — q- |. To simplify the numerical integration, we have given a small width (~ 10 MeV) to the p + q 

nucleon and have used an average for the modulus of p and fixed the (^,^)— plane, avoiding thus to perform two 
nested integrals. The total W~^N — >■ N'p amputated t'^ amplitude is obtained by summing those corresponding to 
the NP, CNP, AP and CAP mechanisms. 

To deduce Eq. (|4ip . we have approximated the ph p— selfenergy (right-hand part of the diagram depicted in Fig. IH]) 

by 

Ticp = {-9c0 + Kh/k") ri\ n^fc) = Fl{k)cJ-^klUx{k) (43) 

with A the charge of the p— meson. Eq. (|43p is obtained by neglecting higher order terms, 0{l'^ /NP), being I — p,qoi 
k. This is consistent with the non-relativistic reduction that leads to the effective potentials in Eqs. (|33p and (|35|) . 
As previously done, we also 

• replace in Eq. (|1T|) . 

By including the non-relativistic Lindhard function for ph + Ah excitations in the denominator, we replace the 
ph excitation of the right-hand in Fig. IHlby a series of RPA excitations through ph and Ah excitations, driven 
by Vt, as depicted in Fig. [TU] (some more details will be given below in the discussion of Eq. ((57)) ). 



multiply t^" by a factor 



t7 ^t^/x{l + (45) 

which allows us to replace the p— exchange interaction in Fig.[S]by the transverse part (Vf) of the effective ph-ph 
potential. 



Finally, and to avoid double counting we must subtract the NP-NP contribution from Eq. (j4T|) . because this term 
was already taken into account in the evaluation of the QE contribution, through the inclusion of a realistic nucleon 
selfenergy, carried out in Ref j2|] . 



This is part of the p— meson propagator, that reads (— g*"^ + k^^k'^ /•m?p)Dp. Actually, only the piece proportional to g^" contributes 



since t^ka = 0. 
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E. The A excitation term 



Like in pion-nuclear, photo-nuclear and electro-nuclear reactions at intermediate energies, the excitation of the 
A(1232) resonance by the weak probe is expected to play a major role. This term is depicted in Fig. (TTJ and its 
contribution to the in medium VF"*"— selfenergy reads 



inw;Ah(9) 



^2cos2 Of 



( 9 



V2%/2 



N 



(27r)4 - + IMaFa Eip ) 



(46) 



with PA — P + Q, M/^ the resonance mass and Fa its width, which can be found, e.g., in Eq. (45) of Ref. (33 |. The 
isospin factor Cn takes the values 1 and -s/S for neutron and proton hole contributions, respectively. Finally, P'^''{pa) 
is the spin 3/2 on-shell projection operator 



P''''{pa) = -{^a + Ma) 



_ 1 _^p1pa. i pir-p^Ar 

^ 3 ' ' 3 M| 3 Ma 



(47) 
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and r°"^{p, q) is the weak A^A transition vertex, that can be found in Eq. (40) of Ref. [321 ■ The contribution to the 
hadron tensor from the selfenergy of Eq. (j46l) is^^ 

Some remarks are here in order. For instance, we must be careful to avoid double counting. Indeed, the contribution 
to the hadron tensor of the Ah excitation term arises from the imaginary part of the A propagator, and in particular 
from the A— width. One of the terms implicit in Eq. (|12p . the one where one picks up the A excitation term both 
in and in jj^* (depicted diagrammatically in Fig. [T2la')'). gives precisely the same contribution plus some medium 
corrections that take into account the Pauli blocking effects. Thus, if we would naively add it to the hadronic tensor, 
the contribution of Eq. (pS]) would be counted twice. Indeed, the term of Fig. fT^a) can be cast in the form of the 



Note that the tensor A'^ can be spUt as in Eq. I I20II . 
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diagram of Fig. [11] but with a A— selfenergy insertion constructed from a pion loop. When the pion is put on the 
mass sheh to build the hadron tensor, we obtain the A— width and thus quahtatively the equivalence is shown^^. In 
a similar way, the diagram of Fig. [T2l b) is one of the terms implicit in the diagram of Fig. |4] that produces a 2p2h 
excitation. 

However, given the iinportance of the A— pole contribution and since the A properties are strongly modified inside 
the nuclear medium a more careful treatment of the A mechanisms is advisable. This implies 

some additional nuclear corrections to Eq. to include the full effect of the self-energy of the A in the medium 
'^a{p{^)) in a systematic manner. In addition, these corrections provide genuinely new contributions to the hadronic 
tensor (e.g. 3p3h mechanisms). Here, we follow the same approach as in Ref. [l|, which is based on Refs. poll4lll48|. 
In the nuclear medium the resonance self-energy is modified because of several effects such as Pauli blocking of the 
final nucleon and absorption processes: AA^ — > NN, AN — > NNtt, or ANN NNN. This is done using a ph-ph 
interaction that includes, besides pion exchange, SRC, a transverse channel driven by p— exchange (see Eq. (|56p ) and 
a RPA-re-summation. 

Following this approach, in the A— propagator, we approximate 



pI-MI+iMaTa V^ + Ma Vi-MA + irA/2 ^ ^ 

with s ~ p\. In the particle propagator of the right hand side of the above equation, we make the substitution: 
rA/2 r^^"i72 - ImSA and take ImSA(/5(r )) and r^^"iV2 as follows. First, the Pauli blockingi^ of the TriV 
decay reduces the Fa free width to F^'^^'', which can be found in Eq. (15) of Ref. Next, the imaginary part 



^'^ More quantitatively, to evaluate this term we can start from Eq. l|lH l and perform the d^fcir integration in order to show explicitly the 

' ^2 [ JuP G'fe Pn) ^ 



Itt A— selfenergy S"''. We get 



in'"' 

W;lphl7r 



N •' 



-iSP° 

A;7r 



(27r)4 |(p + 9)2_M2 -HiMAFAP 

X iTr((;} + M)7«(r-'')t70p„^(-iSi^JPi,r-) (49) 

(pa) = -f:SiVA\' /■ ^£,^(fc^)fc^fe^G'(p + g„fc^.^)(jj + ^_^^^ + j\j) (50) 
\ J J (277)-* 

where /*jya ~ 2.14 is the ttNA coupling, and for simplicity we have evaluated the A— selfenergy for a symmetric Fermi sea. In principle, 
the A— selfenergy which is a matrix in the Dirac space and a Lorentz tensor can be expressed in terms of a linear combination of 

the five orthogonal spin projection operators introduced in Eq. (10) of Ref. [soll . The coefficients, Ai of such linear combination will be 
matrices in the Dirac space and Lorentz scalars. We will enormously simplify the discussion here neglecting A— offshell effects. Within 
this approximation, the spin 3/2 projector of Ref. |50ll reduces to that used here (— P'"'/2Ma) to construct the A— propagator which 
satisfies Pi^^P^s = -2M^P^ and ^aP'^" = P'^''i>£^ = M^^Pf^". Since in Eq. iggjl, E^*^ always appears contracted with two projector 
operators, one realizes that using the orthogonality properties, only the scalar quantity Sa-tt, defined as 

S'i;. = -^A;^P^' + ^ (51) 

will contribute in Eq. I I49I I. which now gets simplified to 



in'"' 



Thus, we recover Eq. 1 148 I I from Eq. I I52I I when we replace G(p;pjv) by i-g^^njv(p)5(p'^ — E{p)) to get rid of the vacuum contribution, 
and consider the imaginary part of the A— selfenergy in Eq. II52I I. This is because to obtain the hadron tensor it appears always ImT,^.^, 
thanks to the symmetry properties of the tensor , and by noting that 

ImEA^. = ^, (53) 

up to density corrections which will account for Pauli blocking. This latter relation follows from the Dyson equation for the A propagator, 
within the on-shell approximation we are using, 

^^A-.TT 72 ^ - 72 »^2 , .»^2^ 1.54) 



In the diagram (a) of Fie. 1121 appears the factor •npf{p)(l — rij^i {p + q — fcjr)), see Eq. II12I I 
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Vind Vi, Vt 

vwv = z z z z : + 




+ ... 



FIG. 13: Diagrammatic representation of the induced interaction. 

of the sclfenergy in Eq. ([55]) accounts for the diagrams depicted in Fig. [T^ where the double dashed hne stands for 
the effective spin-isospin interaction, while the wavy line accounts for the induced interaction. The effective spin- 
isospin interaction is originated by tt and p exchange in the presence of short range correlations. It is obtained by 

substituting o m m 

■'"'^'^ > kfk]k^D^{K) ^ kfk]Vi{k^) + {5,, - k^k])Vtik^) (56) 

with fcf = k'^ /\kTr\. The induced interaction accounts for the series of diagrams depicted in Fig.[T31 There is an RPA 
sum through particle-hole and Ah excitations and it is readily obtained as 



where {/{k^^) is the non-relativistic Lindhard function for ph + Ah excitations^^ (see Eq. (PT|) ). The evaluation of Ea 
is done in Ref. [48]. The imaginary part of Ea can be parametrized as 



, /3 / \ 7 



-ImEA(p(r)) = CQ(^-^j ^^""AjJ ^^^^ 

where the different coefficients are given^° in Eq.(4.5) and Table 2 of Ref. [i^. The separation of terms in Eq. ([55]) is 
useful because the term Cq comes from the AiV — >■ NNtt process (diagrams (c) and (d) of Fig. [T2l when the lines cut 
by the dotted line are placed on shell, and hence the term is related to the {W*,tt) channel), while Ca2,Ca3 come 
from diagrams (b) and (e) and are related to two {W*NN NN) and three (W*NNN NNN) body absorption. 
Hence, the separation in this formula allows us to separate the final cross section into different channels. 

To avoid double counting, we must subtract the contribution of the AP-AP diagram of Fig. I12f a) from Eq. (|19p . 
already taken into account through the r^**"'' piece of the A selfenergy. We must also subtract the contribution of 
the AP-AP diagram of Fig. [T^ b') from Eqs. (pS)) and (|^T|) . because these terms were already taken into account in 
the evaluation of the Ca^ contribution to the A— selfenergy [isj . 



The different coupUngs for A'^ and A are incorporated in f/jv a^nd U/^ and then the same interaction strengths V; and Vt are used for ph 
and Ah excitations 38| . 

The parameterizations are given as a function of the kinetic energy in the laboratory system of a pion that would excite a A with the 
corresponding invariant mass and are valid in the range 85 MeV < < 315 MeV. Below 85 MeV the contributions from Cq and Ca^ 
are rather small and are taken from Ref. [ilj . where the model was extended to low energies. The term with shows a very mild 
energy dependence and we still use the parametrization from Ref. even at low energies. For T„ above 315 MeV we have kept these 
self-energy terms constant and equal to their values at the bound. The uncertainties in these pieces are not very relevant there because 
the A — >■ Ntt decay becomes very large and absolutely dominant. 
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FIG. 14: Irreducible pieces in the Ah channel from the Ah interaction. 



To end this subsection, we would like to devote a few words to the real part of the A— selfenergy and the RPA sum 
of Ah excitations shown in Fig. HM Both of them produce effects on the nuclear response to the weak probe that 
partially cancel out. In Ref. [4^, the dispersive contributions to RcSa associated to the diagrams that gave rise to 
ImEA were also computed. There, it was found ReE^^ ~ —50p/po [MeV] at = 50 MeV and a smooth dependence 
on the pion energy. In principle, ReE^-* could be taken into account by substituting Ma — > Ma + ReE^'' in the 
particle propagator of the right hand side of Eq. ([55)1 . On the other hand, it is easy to realize that the RPA sum 
of Ah excitations, shown in Fig. [T3] can be cast as a contribution to the real part of the A— selfenergy [ij. Actually, 
the latter depends on the particular component of the hadron tensor W^'^ which is being evaluated . Thus, for 
instance, the RPA series depicted in Fig. [TUcan be taken into account, when computing VF^h ^aL (transverse 
components to the direction of the M^— boson) by replacing ReE^^ by ReE^^ + ApVt/Q. This latter sum, in good 
approximation, is positive for the whole range of energies studied here. This was the situation for the inclusive (e, e') 
nuclear reaction studied in Ref. fl| , since there, the excitation of the A resonance by the virtual photon selected the 
transverse mode of the RPA series (see discussion of Eq. (44) in Ref. [ij). However, when the longitudinal component 
^Ah is evaluated, the longitudinal part, Vi, of the effective spin-isospin interaction is selected and now, this RPA sum 
is taken into account by substituting^^ ReE^^ by ReE^^ + ApVi/9, which shows a more pronounced dependence 
than the combination that appeared in the RPA renormalization of transverse components of the hadronic tensor. 
Indeed, it turns out that the ReE^^' + ApVi /9 combination does not have a well defined sign for the whole kinematical 
range of energies studied in this work. Setting to Ma the position of the pole of the A propagator, or changing it by 
adding or subtracting to Ma about 30 MeV, as it could be inferred from the typical values that ReE^-* + 4/9V;(t)/9 
takes for the relevant kinematics to this work, leads to trivial shifts in the position of the A— peak, moderately changes 
of the strength (around 20 %) at the maximum and very tiny changes of the q*^— differential shape. Of course, all 
these effects induced by the RPA-re-summation might be properly taken into account, as it was done for the case 
of the QE-region^^ in Ref. but they, in conjunction with ReE^'', would induce changes smaller than both, the 
precision in the current experimental determination of cross sections, and the uncertainties due to our lack of a precise 
knowledge of the axial nucleon-to- A transition form factor [35| . For simplicity, in this work we will not renormalize 
the real part of the position of the A— peak, which eventually could be studied in the future when more accurate 
measurements become available. 



Note that, in the studies of neutrino induced pion coherent production in nuclei carried out in Refs. l52l |53|| . the replacement ReS^^ 

by ReE^^ + 4pg' /9 is employed to account for the corresponding RPA re-summations. The Landau Migdal parameter g' used there is 
part of Vi, which in addition also includes explicit pion-exchange (see Eq. II33I I) [40l. |41|. 14^1 . This latter contribution was not considered 
in the works of Refs. |52|454|| . because there the distortion of the pion, by using an outgoing solution of the Klein-Gordon equation with 
the optical pion-nucleus potential derived in Ref. [4ll , was implemented, and it accounts for the RPA-renormalization induced by the 
Ah— Ah pion-exchange interaction. 

■^^ For QE kinematics, taking into account properly the RPA effects is much more important 2] than in the A— region, since the cancellation 
of their effects with the difference between particle and hole sclfcncrgies is much less effective. 




FIG. 15: Muon neutrino (left) and antineutrino (right) CC differential cross section ^^'^'^ in oxygen, at 30 degrees of scattering 
angle and with an incident neutrino energy of 750 MeV, plotted against the transferred energy to the nucleus. Different 
contributions are displayed, standing the solid lines for our full model results. Besides in the left panel, we also show results 
(blue dash-dotted line) from Ref. and obtained within the GiBUU framework. 



F. CC antineutrino induced reactions 



The cross section for the antineutrino induced nuclear reaction 

Az ^l+ik')+X (59) 

is easily obtained from the expressions given in the previous subsections, by changing the sign of the antisymmetric 
part of the lepton tensor (L") and using the W^N — > iV'Tr"^ amplitudes of Ref. [32], instead of those involving the 
gauge boson. Note that the pion production off the nucleon amplitudes give rise, directly on indirectly, to all 
contributions considered here, except the QE ones. We take the v—QE cross sections from Ref. [1]. 



III. RESULTS 



We will mainly focus here in the dip and A— peak regions, since the QE contribution was discussed at length in 
Ref. 0. In Fig.[T5l we show results for both muon neutrino (left) and antineutrino (right) induced CC differential 
cross sections at 30 degrees as a function of the energy transferred to the nucleus ("'^^O). The incoming neutrino (anti 
neutrino) energy is 750 MeV. We clearly observe both the A(1232) and the QE peaks; for this scattering angle, the QE 
contribution turns out to be significantly larger than that of the A resonance. We split the full contribution into the 
QE and non QE (A+lplhl7r+ 2p2h) parts. General features are the same for both neutrino and antineutrino induced 
cross sections, and the main difference is an homogeneous reduction in the size of the differential cross section. For 
comparison, in the left panel (blue dashed-dotted line) we also display some results from Ref. [l^, obtained within 
the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) framework, which takes into account various nuclear effects: the 
local density approximation for the nuclear ground state, mean-field potentials, and in-medium spectral functions, but 
does not include those due to RPA correlations. We note first, some discrepancies between these results and ours in 
the QE region, which origin can be traced back to the implementation of RPA corrections in our scheme [3]. Indeed, 
the found differences (small shift in the position and reduction in size, about 25%, of the QE peak) are qualitatively 
identical to those existing between data and GiBUU predictions for the case of inclusive electron cross section for 
a similar kinematics (incident electron energy of 700 MeV and scattering angle of 32 degrees) showed in the upper 
panel of Fig. 9 of Ref. [l^. On the other hand, in this latter figure can be also appreciated the differences with 
the GiBUU model in the description of the dip region. Indeed, we see in Fig. [15] that in the dip region, our model 
predicts larger cross sections than those obtained within the GiBUU scheme. This is due to the 2p2h mechanisms of 
Figs. [31 and [Hi included in our model. Actually, these contributions make also our cross section at the A— peak larger 
than the one predicted in Ref. [T^, even though we use a value of C^(0) smaller than that used in Ref. (1 vs 
1.2). For larger scattering angles, the dip-region cross section becomes relatively much more important, and thus the 




FIG. 16: Muon neutrino (left) and antineutrino (right) CC differential cross section ^^^^ in oxygen, at 60 degrees of scattering 
angle and with an incident neutrino energy of 750 MeV, plotted against the transferred energy to the nucleus. The solid lines 
stand for our full model results. 



inclusion of the 2p2h contributions turns out to be of larger relevance. This is clearly appreciated in Fig. where 
we show results at 60 degrees. In this figure, besides the separation between QE and non QE contributions to the 
differential cross section, the 2p2h part^^ of this latter contribution is shown (orange double dash-dotted curves). The 
blue dash-dotted lines stand in this figure for the results obtained from only the Ah weak-nuclear excitation term 
of Fig. [T^a). neglecting Pauli blocking effects affecting the in medium resonance width. We see how the systematic 
many body VF— absorption modes and the in medium effects considered here change drastically the nuclear response 
fimction in the A— peak, as happened in the QE region as well 0. 

The 2p2h cross section accounts for events were the gauge boson is absorbed by a pair of nucleons, in contrast to 
QE events for which it is absorbed by one nucleon, and furthermore no pions are being produced in this first step. 
Up to re-scattering processes which could eventually produce secondary pions, 2p2h events will give rise to only one 
muon to be detected. Thus, they could be experimentally misidentified as QE events. Yet, IplhlTr events, in which 
the resulting pion from the W absorption is subsequently absorbed and does not come off the nucleus, could be also 
misinterpreted as QE events, if only leptons are being detected. A correct identification of CCQE events, which is the 
signal channel in oscillation experiments, is relevant for neutrino energy reconstruction and thus for the oscillation 
result. By looking at the 2p2h contribution in Fig. [161 we see that at least about 15% of the quasielastic cross section 
might be misidentified in present-day experiments and need to be corrected for by means of event generators. As 
mentioned above, Iplhlvr mechanisms followed by the absorption of the resulting pion, will even make worse the 
situation 'ssj. 

In Fig. [iTl we show CC differential cross sections in carbon for an incident energy of 1 GeV. We observe that 
the 2p2h contribution is sizeable for both, neutrino and antineutrino induced reactions, and that it shows a less 
pronounced dependence than the QE or the A -I- Iplhln components of the total result. On the other hand, the 
antineutrino distribution is much narrower than the neutrino one. Neglecting lepton mass effects, both distributions 
should be equal at q^ = 0, and since the antineutrino cross sections are smaller than the neutrino ones, is reasonable 
to expect the v distributions to be wider than the D ones. 

The MiniBooNE collaboration has measured Q the muon neutrino's CCQE cross section on ^^C. The flux-unfolded 
results as a function of the neutrino energy are depicted in the left panel of Fig. I18[ together with different predictions 
from the scheme presented here. The first observation is that our QE curve misses the data-points, being ourpredicted 
QE cross section significantly smaller than those reported by the MiniBooNE collaboration. Actually in Q, and to 
achieve a reasonable description of the data, an unexpectedly high effective axial mass (entering in the axial- 

vector WNN form- factor) of 1.35 GeV had to be used in the relativistic FG model implemented in the NUANCE event 
generator employed by the MiniBooNE collaboration. This value of Ma is significantly larger than the world average 



A small three body absorption (3p3h) contribution, indueed by Ca^ in the A— selfenergy of Eq. I I58I I is also included under the label 
2p2h in Fig. 1161 and in what follows. 




FIG. 17: Muon neutrino (left) and antineutrino (right) CC differential cross section in carbon for an incident neutrino 
energy of 1 GeV [Q^ = ~<1^)- Different contributions are displayed, standing the solid lines for our full model results. 



value Ma — 1-03 GeV. It is interesting to note, however, that in Ref. [Sj] is also pointed out the NOMAD [23| and 
LSND [5^ high energy {Ei, > 4 GeV) CCQE cross sections are better described with the world average value for Ma- 
The situation become even more worrying, after the work of Ref. [2^ . That work finds that a theoretical approach 
based on the impulse approximation and realistic spectral functions, successfully applied to QE electron scattering, 
fails to reproduce the CCQE neutrino-nucleus cross section, unless the value of the nucleon axial mass resulting from 
deuteron measurements is significantly increased. In addition, they also rule out the possibility, advocated in Ref. 0], 
of interpreting the large Ma resulting from the MiniBooNE analysis as an effective axial mass, modified by nuclear 
effects beyond the FG model [13] • Actually, in [22], it is suggested that the many body techniques successfully applied 
in QE electron-nucleus scattering are not able to explain neutrino induced cross sections and it is argued that the 
development of a new paradigm, suitable for application to processes in which the lepton kinematics is not fully 
determined, will be required. 

Our results do not support this last statement /interpretation, and we rather agree with the picture that emerges 
from the works of M. Martini et al. [27l.[28|. These latter works, in our opinion, constituted a significant step forward 
to clarify the situation. As mentioned above, in the MiniBooNE analysis, ejected nucleons are not detected and the 
QE cross section is defined as the one for processes in which only a muon is detected in the final state. The MiniBooNE 
analysis of the data corrects (through a Monte Carlo estimate) for events, where in the neutrino interaction a pion is 
produced, but it escapes detection because it is reabsorbed in the nucleus, leading to multinucleon emission. However, 
in [23, [15 it is pointed out that 2p2h or 3p3h mechanisms are susceptible to produce an apparent increase in the 
"QE" cross section, since those events will give rise to only one muon to be detected, and the MiniBooNE analysis 
does not correct for them. Within the scheme followed in Ref. j2^], the occurrence of 2p2h final states is described by 
the continuum part of the spectral function, arising from nucleon-nucleon correlations, and there, this contribution 
is found to be quite small (less than 10% of the integrated spectrum). This is not surprising, since our QE results 
(dashed line) in the left panel of FigllS] contain also this contribution^"*, and as we mentioned, we underestimate the 
data. However, the 2p2h contribution considered in [2^ is far from being complete^^ and it corresponds only to the 
many body diagram depicted in Fig. [5] Here, we compute all the contributions contained in the generic diagrams of 
Figs, m [5] and IHl as it was previously done in Ref. [1| for electron scattering, obtained from a realistic model for the 
weak pion production off the nucleon. When these latter contributions are added to the QE prediction of Ref. [§], we 
obtain the solid green line in the left plot of FigUHlin a better agreement with the MiniBooNE data. 

As commented before, these multinucleon knockout events are likely part of the CCQE" cross section measured 



The CCQE cross sections calculated in Ref. were obtained using both particle and hole dressed propagators, determined from a 
realistic in medium nucleon selfenergy |49|| , and thus account for the spectral function effects considered in 22] . 

In fact, the spectral function model taken as a paradigm in the discussion of Ref. 22], though successful to account for the QE electron- 
nucleus scattering, at intermediate energies, badly fails to describe both the dip and the A— regions, as can be appreciated for instance 
in Figs. 5-8 of Ref. Q or in the Fig. 1 of Ref. [22l ]. This is because the lack of a proper model to account for the absorption of the 
virtual photon for two or three nucleons in that model. 
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FIG. 18: Left: Flux-unfolded MiniBooNE z/^ CCQE cross section per neutron as a function of neutrino energy (data-points) 
from Ref. 0], together with different theoretical predictions from this work. Right: Different theoretical calculations for 
antineutrino cross sections per proton off ^^C as a function of the antineutrino energy. For comparison, in both plots, we also 
show the results (blue dash-dotted line) of Martini and collaborators taken from Ref. [2^ . Bands accounting for the theoretical 
errors affecting our results are displayed in both panels. 



by MiniBooNE, and that naturally explains the failure of the scheme of Refs. [sl. [2^. 

Coming back to the left plot of FigfTSl there we also display the band of theoretical uncertainties affecting our results. 
To estimate this band, we have summed in quadratures a 15% relative error in our results with the error induced 
by the uncertainties on the parametrization of the C^{q'^) form factor used here (set IV in Table I of Ref. ^35']). As 
discussed in Ref. [i^ for the CCQE case, 15% is a conservative ansatz to account for the errors, in total and differential 
inclusive cross sections, induced by the uncertainties affecting the nuclear corrections included in our model. Once, 
our theoretical uncertainties are taken into account, we find a reasonable agreement with the MiniBooNE data. We 
would like to stress that we have not fitted here any parameter, and that we have just extended our previous work 
on electron-scattering of Ref. [2| to the study of CCQE cross sections. 

In Fig. [181 we also show the results of Martini and collaborators (blue dash-dotted line), taken from the QE+np-nh 
RPA curves of Fig. 5 of Ref. [1^, which nicely fall within our band of theoretical predictions. Details of the model 
used by M. Martini and collaborators can be found in Ref. [2'^. The evaluation of the nuclear response induced by 
these 2p2h mechanisms carried out in Ref. [2J7|| is approximated, as explained there. The contributions in [27l] that can 
be cast as a A— selfenergy diagram should be quite similar to those obtained here in Subsect. Ill El since in both cases 
the results of Ref. [48] for the A— selfenergy are used. However, some other contributions included here are, either 
not considered or not properly taken into account in [27\ . For example, we believe that none of diagrams of Fig. [5] 
or those in Fig. 2] involving the CT, PP and PF vertices of Fig. [5] have been considered in the work of Martini and 
collaborators. Moreover, the NP-NP, CNP-CNP, NP-CNP, NP-AP, NP-CAP, CNP-AP, CNP-CAP 
and AP — CAP diagrams implicit in Fig. 21 are not directly evaluated in [l^l, but instead, an indirect estimate is given 
for them by relating their contribution to some absorptive part of the p— wave pion-nucleus optical potential. Given 
all this, we find remarkable the agreement exhibited in Fig. 1181 between our results and those previously published in 
Refs. (23,113. 

On the other hand, we see that in our calculation the relative contribution of the 2p2h mechanisms with respect to 
the QE cross section, is quite similar for both neutrino and antineutrino induced processes. Thus, in what respect to 
this issue, our results do not support the claims of Ref. [l^ on a minor role of the 2p2h mechanisms in the antineutrino 
mode. 

We should mention that the MiniBooNE collaboration has also published the flux-integrated CCQE" double 
differential cross section cP a / dE^d cos 9 fj, in bins of muon energy E^ and cosine of the muon scattering angle with 
respect to the incoming neutrino direction. We must refrain to compare with these valuable data. The reason is 
that the MiniBooNE flux remains sizeable up to neutrino energies too high to make meaningful the predictions of 
the model presented here. Indeed, neutrino energies of 1 or 1.2 GeV at most, is the clear upper limit of validity of 
our predictions. The fraction of the MiniBooNE flux above 1.2 GeV is still larger than 17%, and this together with 
the fact that the cross sections grows with the energy has prevented us to do the comparison. We are working in 
extending our model to higher energies, but this is far from being a trivial task. New channels (2p2hl7r or two pion 
production, . . . ) and higher resonances must be incorporated into the model, besides of smoothly getting rid of the 




FIG. 19: Left: Data points stand for the SciBooNE neutrino CC inclusive interaction cross section per nucleon f6|. We also 
show our QE and full model results, and in this latter case the theoretical uncertainty band is also displayed. At 1 GeV, we 
depict the full model cross sections obtained when the GTR value of 1.2 for C^{0) is used instead of 1 (violet triangle), and 
when some 2p2hl7r contributions (blue empty circle), neglected in the present work, are taken into account (see text for some 
more details). We also give at 1.1 GeV the total cross section obtained with C^(0) — 1.2. Right: QE and full model predicted 
antineutrino CC inclusive cross section per nucleon, as a function of the antineutrino energy. 



RPA and other nuclear effects, which importance diminishes when the transferred energy increases. 

To end this section in Fig. [19] we show total and QE inclusive cross sections for both neutrino and antineutrino 
modes. In the neutrino case, we compare our results with the recent data published by the SciBooNE collaboration. 
We display SciBooNE data-sets based on NEUT and NUANCE Monte Carlo event generators. We find a reasonable 
description, taking into account experimental and theoretical uncertainties, up to neutrino energies around 1 GeV. 
At larger energies, we underestimate the cross section, as anticipated above. For instance, we see how some 2p2hl7r 
contributions neglected in our model, become relatively important a.t = 1 GeV. More specifically, the empty circle 
is obtained when the A— resonance contribution to the imaginary part of U is kept in the evaluation of the imaginary 
part of the it— and p— selfenergies in Eqs. ((29|) and (|44p . There are some other W'^NN — ^ NNtt mechanisms which 
should be taken into account, as well as the contribution of higher resonances ^]. Though small, also kaon [59| . 
hyperon [6^ and two pion [61] production channels should be considered to end up with a robust theoretical model 
above 1 GeV. 



IV. CONCLUSIONS 



We have developed a model for the study of weak CC induced nuclear reactions at intermediate energies of interest 
for current and future neutrino oscillation experiments. This model is an extension of the work of Ref. that analyzed 
the QE contribution to the inclusive neutrino scattering on nuclei. The model is based on a systematic many body 
expansion of the gauge boson absorption modes that includes one, two and even three body mechanisms, as well as 
the excitation of A isobars. The whole scheme has no free parameters, besides those previously adjusted to the weak 
pion production off the nucleon cross sections in the deuteron, since all nuclear effects were set up in previous studies 
of photon, electron and pion interactions with nuclei. 

We have discussed at length the recent CCQE MiniBooNE cross section data. To understand these measurements, 
it turns out to be essential the consideration of mechanisms where the M^— boson is absorbed by two or more nucleons 
without producing pions, as first suggested by M. Martini and collaborators [27j|. Our evaluation of these pionless 
multinucleon emission contributions to the cross section is fully microscopical and it contains terms, which were either 
not considered or only approximately taken into account in [27j . We end up with a reasonable description of the 
neutrino CCQE" MiniBooNE and total inclusive SciBooNE cross section data up to neutrino energies of around 1 
GeV. 

Our results do not support the incompatibility among neutrino and electron-nucleus inclusive data claimed in 
[2^ , since our neutrino model is just a natural extension of that developed in Refs. [ij and [s^l to study electron- 
and photo-nuclear inclusive reactions. Indeed, we believe that the origin of the problem can be traced back to the 
difficulties of the spectral function model advocated in [22| to properly describe the dip and A— regions for electron 



24 



scattering, together with the mismatch existing in the definition of the quasielastic contribution between the theory 
and the experimental neutrino communities. 

The recent CC MiniBooNE and SciBooNE inchisive data sets provide very valuable information to distinguish 
among different models. This will definitely help to unravel the details about the modification of the CC weak cur- 
rent properties inside of the nucleus, and will set up the basis to construct a robust theoretical framework where 
all electroweak nuclear reactions at intermediate energies could be studied. This is in turn of special relevance to 
better understand the systematic errors affecting present (MiniBooNE & T2K) and coming neutrino oscillation ex- 
periments involving neutrinos with energies below 1 GeV. Future antineutrino data, similar to the CCQE MiniBooNE 
measurements of total and differential neutrino cross sections, will further constraint any theory. 

We think the microscopical model presented here, which extends that of Ref. 2] beyond the QE region, constitutes 
a first step towards this goal. The model should be extended still to higher energies, which would make possible the 
comparison of its predictions to the MiniBooNE differential cross section data. This is a non trivial task, but it will 
allow for a better knowledge of the axial currents both for hadrons and nuclei. 
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